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Abstract. Recent studies of the substellar pop ulation in the Taurus cloud have revealed a deficit of brown dwarfs 
comp ared to the Trapezium cluster population fericefio et aljlT998t lLuhmanll20o"oT : iLuhman et alJl2003at ILuhmanl 
2004). However, these works have concentrated on the highest stellar density regions of the Taurus cloud. We 
have performed a large scale optical survey of this region, covering a total area of ~ 28deg 2 , and encompassing 
the densest parts of the cloud as well as their surroundings, down to a mass detection limits of 15 Mj. In 
this paper, we present the optical spectroscopic follow-up observations of 97 photometrically selected potential 
new low-mass Taurus members, of which 27 are strong late-M spectral type (SpT > M4V) candidates. Our 
spectroscopic survey is 87 % complete down to i! —20 for spectral types later than M4V, which corresponds to 
a mass completeness limit of 30 Mj for ages < 10 Myr and Av < 4. We derive spectral types, visual absorption 
and luminosity class estimates and discuss our criteria to assess Taurus membership. These observations reveal 
5 new VLM Taurus members and 12 new BDs. Two of the new VLM sources and four of the new substellar 
members exhibit accretion/outflow signatures similar to the higher mass classical T Tauri stars. From levels of 
Ha emission we derive a fraction of accreting sources of 42 % in the substellar Taurus population. Combining 
our observations with previously published results, we derive an updated substellar to stellar ratio in Taurus of 
IZss = 0.23 ± 0.05. This ratio now appears consistent with the value previously derived in the Trapezium cluster 
under similar assumptions of 0.26 ±0.04. We find indication that the relative numbers of BDs with respect to stars 
is decreased by a factor 2 in the central regions of the aggregates with respect to the more distributed population. 
Our findings are best explained in the context of the embryo-ejection model where brown dwarfs originate from 
dynamical interactions in small N unstable multiple systems. 

Key words. Stars: low-mass, brown dwarfs - Stars: late-type - Stars: luminosity function, mass function - Stars: 
pre-main sequence 



1. Introduction 

Recent works bv IBriceno et all (ll99Sl ): ILuhmanl feOOdk 
IBriceno et all l)2002|) ; ILuhman et alJ (|2003a|k ILuhmanl 
l)2004|h have revealed a factor of 1.4 to 1.8 deficit of 
Brown Dwarfs (BDs) in the Taurus cloud compared to the 
Trapezium cluster, possibly indicating that sub-stellar ob- 
ject formation depends on the environment. If confirmed, 
the lower abundance of BDs in Taurus could have strong 
implications on substellar Initial Mass Function (IMF) 
models. However, as these previous studies were concen- 
trated on high stellar density regions, the majority of the 
volume occupied by the molecular clouds was left unex- 
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plored. There is a possibility that Taurus BDs are not 
clustered in the same regions as most of the T Tauri stars, 
or that they have scattered away from their birth sites as 
proposed bv lReipurth fe Clarke! l|200l|) . If BDs are ejected 
during their formation, a significant fraction of the substel- 
lar content of the central parts of the cloud may have been 
missed. Indeed, with ejection velocities > 1 km.s" 1 , as pre- 
dicted by some of the ejection m odels ijReipurth fc Clarke) 
l200lUKroupa k Bouvierl2003a|h BDs could have travelled 
as far as 1° from their birth sites at the Taurus distance 
in 2 - 3 Myr. 

The recent availability of large visible and infrared 
cameras makes it possible to quickly survey large por- 
tions of the sky encompassing tens of square degrees down 
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to the sensitivity needed to study the substellar domain. 
In order to check for the existence of a widespread BD 
population in Taurus, we have used the Canada-France- 
Hawaii telescope to perform a large optical survey covering 
a total of 28deg 2 down to completeness limits of z'=21.8, 
z'=20.9, corresponding to a mass detection limit of 15 Mj 
for an age < 5 Myr and Ay < 5, according to the pre- 
mam sequence DUSTY model of lChabrier etaTI $20()(t\ . In 
this paper, we report on medium resolution optical spec- 
troscopic follow-up observations of 97 photometrically se- 
lected Taurus candidate members, of which 27 are strong 
late-M (SpT > MW) candidates. This work is the con- 
tinuation and com pletion of the prelim inary spectroscopic 
study published in lMartin et alJ l)200l|) . where initial can- 
didates were selected from a 3.6 deg 2 photometric survey 
conducted with the CFHT12k camera, which is a subset 
of the full 28 deg 2 CFHT Taurus survey described here. 

In section 2, we first briefly describe the CFHT pho- 
tometric survey and candidate selection procedures, then 
the spectroscopic observations presented here. We derive 
spectral types, reddenings and luminosity class estimates 
from both a spectral fitting procedure and the analysis of 
colors and Na I equivalent widths in section 3. We assess 
Taurus membership in section 4 and discuss the properties 
of the new 5 VLM members and 12 BDs identified in this 
study. In particular, we investigate their accretion/outflow 
signatures and spatial distribution and update the stellar 
to brown dwarf ratio in Taurus. Finally, we discuss the 
implications of our results on the models for substellar 
object formation in section 5 and summarize our findings 
in section 6. 



2. Observations 

2.1. The CFHT optical photometric survey 

The overall Taurus survey has been performed in 4 succes- 
sive periods on the Canada-France-Hawaii telescope with 
the CFH12k and MEGACAM large-scale optical cameras 
(see Tabled for a detailed journal). The technical char- 
acteristics of the CFH12k and MEG ACA M cameras are 
presen ted in ICuillandre et al.1 l)2000|) and IBoulade et alJ 
i|2003|) respectively The first part of the CFH12k survey 
(centered on the densest part of the Taurus cloud, from 
1999 to 2001), has been obtained as part of CFHT di- 
rector's discretionary time, while the remaining larger set 
of data has been obtained, in service mode, as part of a 
larger key program devoted to the study of young clusters. 
A detailed description of the optical photometric survey is 
deferred to a forthcoming publication (Guieu et al 2005 in 
preparation). We summarize here its main characteristics. 

Data reduction , performed at CFHT using e lements of 
the Elixir system l|Magnier fc Cuillandrell2004|) . included 
bias and dark subtraction, flat-fielding, fringing correc- 
tion, bad pixels removal and individual frame combina- 
tion. Point source detection was performed on the com- 
bined I + z' images. For the CFH12k data, PSF fitting 
photometry was extracted with the PSFex routine from 



the SExtractor program I Berlin fc Arnoutslll996h . while 
aperture photometry was obtained for the MEGACAM 
data with the same program. Photometric catalogs were 
combined, using the following transformation between 
CFH12k (I,Z') and MEGACAM (i',z') photometric sys- 
tems, computed from overlapping fields: 



and 



i' = I + 0.85 * (I - Z') + 0.23 



z' = 0.25 * J + 0.746 *Z' + 0.32 



The total survey yielded more than 10 6 sources de- 
tected down to i'=2A and z'=23. From the turn-over at 
the faint end of the magnitude distribution, we estimate 
the completeness limits of our optical photometric sur- 
vey to be i'=21.8 and z'=20.9. Concerning the brightest 
sources in the survey, the saturation limits are i' = 12.5 
and z' = 12 

We identify Taurus low-mass candidate members with 
the following procedure. We first select optical sources 
with a detection in the near-infrared 2MASS catalogue 
at both J, H and K bands with quality flags A, B or C. 
We then evaluate a rough estimate of the reddening vec- 
tor towards each source by assuming a maximum intrinsic 
(J — H) color of 0.8, which correspo nds to a spectral 
type M9-L0 for a main sequence dwarf ^Kirkpatrick et all 
l2000HDahn et al.ll2002t) . Reddening values estimated this 
way will be conservative lower estimates for earlier spec- 
tral type sources. Candidate Taurus members are then 
selected according to their location in the various color- 
magnitude and color-color diagrams plotted in Figure d 
We dereddened each candidate and required its position in 
the —z') (or I/(R — I) when available) and — J) 
color-magnitude diagrams to be compatible with stars 
younger than 10 Myr , derived from comparison with the 
IChabrier et all §000) pre-main sequence tracks. We also 
require colors compatible with a dwarf M spectral type 
in the (i' — J) /(J — K) diagram. This selection procedure 
yielded 264 Taurus low-mass candidates, of which 227 have 
i' band magnitudes <20 . 

Strong galactic contamination from primarily early M- 
type background giants is expected. In order to properly 
assess Taurus membership, spectroscopic follow-up of the 
photometrically selected candidates is therefore manda- 
tory. We restrict ourselves to the 227 sources with i' band 
magnitudes < 20, in order to allow a proper spectro- 
scopic follow-up in the optical domain. We further use the 
(J—H) I (H—K) diagram to separate early and late M can- 
didates. This procedure also allows to minimize the back- 
ground giant contamination, expected to be very strong at 
the early M spectral types. Among the 227 sources with i' 
band magnitudes <20, 47 sources appear later than M4V 
in the (J - H)/(H - K) (see Figured!). 

Among t hese 47 sources , are in cluded 4 sources already 
published in iMartm et"afl l|200lh (CFHT-Tau 1, CFHT- 
Tau 2, CFHT-Tau 3 and CFHT-Tau 4) and 10 previously 
known Taurus low-mass members (V410 Tau Anon 13, 
ZZ Tau IRS, J04442713+2512164, J04381486+2611399, 
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Fig. 1. Observed (not dereddened) color-magnitude and color-color diagrams used to select low- mass Taurus candi- 
dates. Small black dots in panels a), b) and c) are the 264 sources that pass through the selection process (see text 
for more details). Only the 227 sources with i' < 20. are represented in panel d). The 44 (40 width i' < 20) primary 
photometric candidates observed spectroscopically are displayed in filled red circles. Open blue circles are previously 
kno wn Taurus members . Green squares identify the spectroscopically confirmed new Taurus members from this work 
and iMartm et all ll200l|). In panels a ) and b ), the full black curve shows the location of the 10 Myr isochrone from 
the DUSTY model of lchabrier et alJ l)2000|) at the Taurus distance. In panel c), the black curve shows the straight 
line approximating the dwarf sequence used to select M-type candidates. The arrow indicates a reddening vector of 
Ay = 4 magnitudes. In the (J — H)/(H — K) diagram are also shown the locations of the dwarf (green curve) and 
giant (open red triangles) branches, as well as reddening vectors extending from the M4 and M6 dwarf colors (dashed 
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Instrument 


FOV 


Pixel 


Date 


Area 


Band 


Detection 


Mass limit 




(deg 2 ) 


("/pixel) 




(deg 2 ) 




limit 


{A v = 5, 5Myr) 


CFHT12k 


0.33 


0.21 


1999-2001 


3.6 


R, I, Z' 


23, 22, 21 


15Mj 


CFHT12k 


0.33 


0.21 


2002 


8.8 


I, Z' 


22, 21 


15Mj 


Megacam 


1 


0.19 


2003-2004 


24 


i\ z' 


24, 23 


10M.j 



Table 1. Overview of the CFHT optical survey of the Taurus cloud 



KPNO-Tau 9, KPNO-Tau 7, KPNO-TAU 5, KPNO- 
Tau 3, J04284263+2714039, J04202555+2700355). We 
present in this paper medium resolution (R~ 1000) optical 
spectroscopic observations obtained for 27 of the remain- 
ing 33 unknown late-M (SpT > M4V) can didates includ- 
ing K P NO Tau 1-2-4 and 6 p ublished in IBricefio et alJ 
l)2002|) : iLuhman et"aTl (|2003a() . Thus, combining obser- 
vations presented in this work and previously published 
results, spectroscopic information is available for 87 % 
(41/47) of the Taurus candidates later then M4V iden- 
tified from our CFHT survey. We also included in our 
spectroscopic survey 13 candidates with SpT < MAV, 
which represent 7 % (13/180) of the remaining early M- 
type candidates with i' < 20. A very strong galactic con- 
tamination is expected in this latter sample. We note that 
our spectroscopy survey is 100 % complete in the sub- 
stellar domain: all the 13 sources that appear later than 
M6V in the (J — H)/(H — K) diagram are included in our 
spectroscopic study. 

In addition to the 44 (40 with i' < 20) sources selected 
with the procedure just described, we have also observed 
in the course of the different spectroscopic campaigns, 
53 additional sources selected with less restrictive crite- 
ria. These include in particular sources selected from the 
I/(R — I) diagram, obse rved at the WHT i n September 
2000 but not published in lMartfn et~all l)200l|) . All of these 
additional objects, turn out to be spectroscopically con- 
firmed as non Taurus members (see They are how- 
ever included in this paper and their derived properties 
are summarized in Table |3J 

Due to our selection procedure, the mass-completeness 
limit of our spectroscopic sample is set both by the 2MASS 
completeness limits (J=15.25, H=14.4, K=13.75) and the 
sensitivity limit of our optical spectroscopic observations 
(i' <20). Both are equally restrictive in terms of mass 
completeness (see M4.4JI , 

2.2. Spectroscopy 

Spectroscopic observations have been conducted wit h the 
WHT/ISIS, FORS1/VLT llAnpenzeller et all Il998l) and 
LRIS/Keck instruments over 5 runs between September 
2000 and December 2004. The detailed journal of the spec- 
troscopic observations is presented in Table HJ 

2.2.1. WHT/ISIS observations 

Spectroscopic observations were carried out using 
the Intermediate-dispersion Spectrograph and Imaging 



System (ISIS) at the 4.2 m William Herschel Telescope 
(WHT) in La Palma on 2000 September 2829. The R158R 
grating on ISIS's red arm gave a wavelength coverage from 
640.9 to 936.5 nm. The spectral resolution was 2.5 pixels 
(7.2 ). The data were reduced using standard routines for 
bias subtraction and flat-field correction within the IRAF7 
environment. Wavelength calibration was made using the 
spectrum of an NeAr lamp. Instrumental response was cal- 
ibrated out using spectra of the flux standard Feige 24. A 
total of 29 source s were observed durin g this campaign, 4 
were published in iMartfn et al.l l)200lj) and 6 turn out to 
be already known Taurus members. We publish here the 
remaining 19 sources. 



2.2.2. VLT/FORS observations 

FORS1 1 observations were performed in service mode be- 
tween 2001 and 2003. In both cases, the grism 3001 was 
used, providing a dispersion of 2.59 A/pixel and a wave- 
length coverage of 6000 to 11000 A. The slit width was 
0.7". The standard ESO pipeline reduction procedure has 
been applied. This includes bias correction, flat-fielding 
and wavelength calibration. Individual spectra were ex- 
tracted using the IRAF 2 apextract routine. The spectral 
response was calibrated with observ ations of the standar d 
M9 field BD DENIS-P-J1048-3956 l|Delfosse et al.ll200l|) . 



2.2.3. Keck/LRIS observations 

Two observing runs were conducted with the LRIS /Keck I 
spectrograph in January 2003 and December 2004. The 
spectral dispersion was 1.87 A/pixel and wavelength cov- 
erage 6000-10000 A. The slit width was 1.0" for the first 
run and 1.5" for the second one, due to worse seeing 
conditions. Individual exposures were first bias corrected 
and flat-fielded. Spectra were extracted with the IRAF 
apextract routine, then wavelength calibrated with a ref- 
erence arc lamp exposure. The spectral response is cali- 
brated with observations of the following standards: LTT 
377, Feige 110, BD+28 4211, BD+25 4655, PG 0918+029, 
G193-74, G191B2B, Feige 11. 



1 FOcal Reducer/low dispersion Spectrograph 

2 IRAF distributed by National Optical Observatories 
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id 


Instrument 


Telescope 


Dates 


Wavelength 


Dispersion 


Number of sources 








A 


A/pixel 






1 


FORS1 


VLT-UT3 


2001-10-25 to 2002-02-20 


6000-11000 


2.59 


10 


2 


FORS1 


VLT-UT1 


2002-12-11 to 2003-03-05 


6000-11000 


2.59 


10 


3 


LRIS 


Keck I 


2003-01-08, 2003-01-09 


6000-10000 


1.87 


26 


4 


LRIS 


Keck I 


2004-12-03, 2004-12-04 


6000-10000 


1.87 


32 


5 


ISIS 


WHT 


2000-09-28, 2000-09-29 


6490-9365 


2.88 


19 



Table 2. Journal of spectroscopic observations 




4 6 8 10 7000 7500 8000 8500 

SpT (M) X (A) 

Fig. 2. Left : x 2 map obtained in the Ay-SpT plane for one of our candidate spectrum (2MASS J04390396+2544264). 
In solid lines are also plotted the Ay variation vs SpT derived from the (I-J) and (J-H) colors using the 
iRleke fc Lebofskvl l|l985h extinction law. Dashed lines correspond to 1 a uncertainties. The best fit obtained gives 
an M7.25 spectral type and Ay — 0.4 mag, compatible with the broad band colors. The best fit is obtained for the 
average of a dwarf and a giant template, indicating a luminosity class IV. Right: The resulting best fit ( average of 
a dwarf and giant template of SpT M7.25) is represented with a dashed red curve overlaid on the observed spectrum 
(solid black curve). 



3. Data analysis 

3.1. Spectral-type and reddening determination 
3.1.1. Numerical spectral fitting 

Spectral types can be derived from spectro-photometri c 
indices, such as the ones defined bv iMartm et alJ (^)99). 
However, these indices are affected by the strong extinc- 
tion expected towards Taurus members. From simulations 
of artificially reddened standard spectra, we have com- 
puted that for a visual extinction of 5 mag nitudes, spec- 
tral t ypes estimated from the PC3 indice 1 Martin et alJ 
1999) are wrong by 2 classes for types around M7V. We 
therefore chose to apply a spectral fitting proce dure simi- 
lar to the one presented in lLuhman et alj l)l998[K allowing 
to derive simultaneously spectral type and reddening. We 
built a library of reference dwarf and giant spectra span- 
ning the range between Ml and L0. The dwarfs reference 
library w as built from the c ombin a tion of spectra pub- 
lish ed in iKirkpatrick et~afl l)l999h . IMartm et all l|l999|) 
and iLeeeett et alJ l)2000|) . with 4 individual spectra, on 
average, for each M-spectral subclass. The giant reference 
library, built fro m a sample of 97 y ery bright M-giants, 
was taken from iFluks et al] l|l994h . Half and quarter 
sub-classes were constructed by combination of adjacent 
classes (for instance M6.25 =< 0.75 x A/6 + 0.25 x M7 >). 



For a given spectral type, each candidate spectrum can be 
fitted either by the dwarf template, the giant template or 
the average of both. This procedure therefore allows for 
an estimate of the luminosity class. 

The fit is performed in the wavelength range 7000- 
8500 A, where the signal to noise is maximized and the 
spectral response is best determined. Figure shows the 
X 2 map in the {Ay - spectral type) plane for one of our 
Taurus candidate (2MASS J04390396+2544264) with the 
corresponding best fit solution. In the % 2 map, we have 
also added the Ay versus spectral type relations pre- 
dicted by the (I-J) and (J-H) colors. The calibration rela- 
tions of broad band colors versus spectral type are derived 
from linear fits to the photometr y of M - type dwarf stan- 
dards observed by iLeggett et all <|l996L Il998l l200ltl and 



iDahn et alJ 1 2002fl . In the example shown in Figure El the 
best fit solution found by our numerical fitting procedure 
appears compatible, within the estimated lcr uncertain- 
ties, with the solution predicted by the (I-J) and (J-H) 
colors. From simulations of artificially reddened standard 
star spectra, we estimate typical internal uncertainties for 
our numerical fitting procedure of half a class in spectral 
type and 0.8 magnitudes in Ay. Among our initial spec- 
troscopic sample, 33 sources are classified as giants, 38 
as field dwarfs and 26 as sub-giants (ie best represented 
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by the average of a dwarf template and giant template). 
These 26 sub-giants appear as our primary Taurus pre- 
main sequence candidates. We investigate below in more 
details the reliability of the parameters, in particular red- 
denings and luminosity classes, determined by the spectral 
fitting procedure. 

3.1.2. Reddening estimates from color-color diagrams 

In order to assess in more details the reliability of the 
Ay values derived from our numerical fitting procedure, 
we independently estimate extinction values from the (J- 

H) /(H-K) and (I-J)/(J-K) color-color diagrams. We ex- 
clude from the following analysis the 33 sources classified 
as giants by the spectral fitting procedure. Most of these 
turn out to be early type (SpT < M2III) background gi- 
ants. We are therefore left with a sample of 64 sources 
classified either as M dwarfs or sub-giants. Sources with 
no detectable excess infrared emission in the (J-H) versus 
(H-K) diagram are dereddened back to the locus of dwarfs, 
while those lying to the right of the reddened dwarf se- 
quence in the near-IR color-color diagram are dereddened 
back to th e locus of ac t ively accreting T Tauri stars, as 
defined bv iMever et"aD l|l997h . In the (I-J) versus (J-K) 
diagram, sources are dereddened back to the dwarf se- 
quence. 

We compare in Figure 01 (Left) Ay values derived in 
these two ways. Four sources (individually labelled) show 
large differences in excess of 2 magnitudes. Two of these 
sources also show significant near-infrared excess emission 
in the (J-H) / (H-K) diagram (Figure [SJ , indicating that 
the determination of Ay is likely very uncertain. Except 
for these 4 sources, Ay values derived from both color- 
color diagrams agree remarkably well: the distribution of 
differences peaks at AAy=0 with a dispersion of 0.3 mag- 
nitudes. On the other hand, we find a systematic difference 
between the Ay values obtained from the spectral fitting 
procedure and the ones derived from the color-color dia- 
grams (Figure E3 Right), especially for the sources classi- 
fied as sub-giants: the distribution of differences peaks at 
AAy=l with a dispersion of 0.3 magnitudes. For sources 
classified as dwarfs the distribution of differences peaks at 
A/l\/=0.25 with a dispersion of 0.3 magnitudes. 

The discrepancy observed in Ay determinations for 
the sub-giant population likely arises from an improper 
calibration of the intrinsic broad-band colors for sources 
with surface gravities intermediate between dwarfs and gi- 
ants. Indeed the derivation of reddening values from main- 
sequence colors assumes that a young object has the same 
intrins ic continuum slope as a dwarf. However, iLuhmanl 

I) l999|) noted that intrinsic (I — J) colors of pre-main se- 
quence sources are significantly redder than the ones of 
dwarfs at spectral types later than M5V. Indeed, the spec- 
tra of the mid to late M young stellar objects always rise 
more rapidly beyond 8000A than the spectra of main- 
sequence dwarfs of the same spectral type. The deriva- 
tion of reddening values from dwarf colors involving the 



I band will therefore systematically overestimate Ay for 
young mid to late-M stars. Our analysis seems to indicate 
that this is also true for reddening values derived from 
near-infrared colors. The spectral fitting procedure on the 
other hand takes into account the fact that the surface 
gravity of young stars is intermediate between the ones of 
dwarfs and giants. Although we note that the continuum 
slope of a young star may not be exactly the average of 
a dwarf and giant continuum, we estimate that the Av 
values derived from the spectral fitting procedure are less 
biased and we adopt them in the following. 

3.2. Luminosity class and Sodium equivalent widths 




2000 2500 3000 3500 



Teff (K) 

Fig. 4. Measured EW(Na 1) vs effective temperature for 
our spectroscopic sample. Star symbols represent spectra 
best fitted with the average of a giant and field dwarf 
template (sources classified as luminosity class IV). tri- 
angles show sources best fitted with a field dwarf tem- 
plate; circles correspond to a giant template. Sources with 
EW(Ha) > 8 A are plotted in red. Solid lines show pre- 
dictio ns from the DUSTY models from IChabrier et alJ 
(2000j spanning log 10 (g/go) =3.5 to 6 by step of 0.5 (as la- 
belled); dashed lines correspond to the COND model span- 
ning log 10 (g/go) =2.5 to 6. The 4 sources with discrepant 
EW(Na 1) are indicated with an arrow. 

The Nai doublet at 8183 and 8195 A is sensitive to 
surface gravity effects ijMartm et alJ fl~999) for effective 
temperatures in the range 2500 to 3500 K. We therefore 
compare the luminosity class estimated for our sources 
from the spectral fitting procedure with Na 1 equivalent 
widths measurements. The EW(Nai) is computed between 
8172 A and 8207 A (the doublet is not fully resolved at 
the moderate spectral resolution of our observations). The 
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Fig. 3. Comparison of Ay values derived from the color-color diagrams and the spectral fitting procedure. Left: 
Histogram of differences in Ay values derived from the (J-H)/(H-K) and from the (I-J)/(J-K) diagrams. Sources 
classified as giants have been excluded. Right: Histograms of differences in Ay values derived from the (I- J) / ( J-K) 
colors and from the spectral fitting procedure, in full line for sources classified as sub-giants, in dotted line for sources 
classified as field dwarfs. Individual sources with the largest Ay discrepancies are labelled and discussed in the text. 



local continuum is estimated by a linear fit to the three 
highest points in each of the following spectral intervals: 
[8090:8172 A] and [8207:8280 A]. We estimate typical un- 
certainties on the equivalent widths to be on the order of 
5 %. 

Figure 01 illustrates the excellent agreement between 
the luminosity class estimated from our spectral fitting 
procedure and the Nai equivalent widths. In this figure, 
we plot the measured EW(Nai) versus effective tempera- 
ture T e ff. Spectral types, derived from the spectral fitting 
procedure, are converted in to effectiv e temp eratures using 
the relationships given by iLuhmanl (^99) for pre-main 
sequence stars, dwarfs and giants respectively. Eighteen 
sources classified as early- type (SpT < Ml) giants are 
not included in the graph. Figure 0] shows a clear sepa- 
ration between sources spectroscopically classified as field 
dwarfs and sources identified as sub-giants. All the sources 
classified as field dwarfs show EW(Nal) co mparable to the 
predic tions of the model atmospheres from lChabrier et alJ 
( 2000J) with surface gravities log 10 (g/go) > 4. On the other 
hand, all but 4 of the estimated luminosity class IV sources 
show EW(Nai) consistent with predictions from the model 
atmospheres for low surface gravities (log 10 (g/go) < 4). In 
the following, we therefore adopt the luminosity class es- 
timated from the spectral fitting procedure and discuss 
the nature of these 4 discrepant low surface gravity candi- 
dates (indicated by an arrow in Figure^}. We notice that 
one known Taurus BD (KPNO-Tau 4) shows an unusu- 
ally small Nai equivalent width (2.07 A) for its estimated 
effective temperature (< 2500 K). However, the detailed 
profile of the Nai doublet in this source appears to differ 
significantly from both the model predictions and the cor- 
responding dwarf template line profile, suggesting a core 
emission which could explain its low position in Figure 0] 



4. Results 

4.1. Identification of Taurus members 

We compile in Table El and 0] the reddening, spectral type 
and luminosity class derived from the spectral fitting for 
our 97 spectroscopic sources. We also list measured equiv- 
alent widths in Nai and Ha. We use the luminosity class 
as the primary criterion to identify Taurus pre-main se- 
quence sources. 26 sources are classified as luminosity class 
IV through our spectral fitting procedure. Among these, 
only 4 show EW(Na i) indicative of larger surface gravities 
log 10 (g/g ) > 4, as determined from Figure^ 

The first one, 2MASS J04221675+2654570, has an 
early spectral type (Ml. 25), and falls in the tempera- 
ture regime where the sodium equivalent width is no 
longer sensitive to surface gravity for moderate values 
(log 10 (g/g ) < 5). In addition, this source shows signif- 
icant Ha emission (with an equivalent width of -24.5 A), 
reddening (Ay — 6.6) and near-infrared excess (see Fig. 
IHJ). We therefore classify it as a Taurus member. 

The second is 2MASS J04355760+2253574. The spec- 
trum of this object is best fitted by the average of a M4.75 
dwarf and giant with a reddening of Ay — 2.55 and it dis- 
plays Ha emission with an equivalent width of -21.58 A, 
indicating a possible pre-main sequence nature. Its near- 
infrared colors ((J-H)=1.01 and (H-K)= 0.56) suggest sig- 
nificant reddening and/or near- infrared excess emission. 
However, its measured EW(Nai) of 5.25 A is more rep- 
resentative of field dwarf surface gravities. In addition, 
when placed on the HR diagram at the Taurus distance 
with the above derived parameters, it falls close to the 
100 Myr isochrone, significantly below the Taurus popu- 
lation. The spectral fit obtained with a dwarf template 
gives a spectral type of M5.25 and Ay of 0.45 with a x 2 
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2MASS Name 


Run RA(J2000) DEC(,J2000) Av 


EW(Na) EW(fl„) 


SpT 


Field stars 



1 <"\ Inl ^ in Oi ; Or 1 (hi ■ 

j u4uio 1 4u-t--£b.£biyy * 


4 


04 


01 


C7 
O ( 


409 


26 


2b 


19 


993 


1 .59 


4.20 


. 5 5 


_' 


00 


V 


J 041oOo88+28Uoooy* 


1 


04 


13 


08 


887 


28 


05 


55 


943 


0.00 


8.39 


1.69 


6 


75 


V 




5 


04 


14 


19 


283 


28 


04 


1 1 


813 


0.47 


6.18 


0.58 


3 


50 


V 


J U4141yOl + z ; CH U14* 


1 


04 


14 


19 


519 


27 


57 


01 


45 1 


0.90 


8.60 


-7.90 


5 


75 


V 


TH/l 1 7(11 7Q 1 1001 EflQi 

JU41 t Ul /o+ziozloyo* 


5 


04 


17 


01 


787 


28 


2 1 


59 


375 


8.39 


3.62 


0.00 


1 


00 


III 


iri/i i 7Q7/IC i noiinQn, 
J U41 (oi 40-1-^011^00* 


5 


04 


17 


37 


463 


28 


1 1 


23 


046 


9.82 


2.77 


0.00 


1 


00 


III 


Tn/i i onm r: i noi 7107 , 
J 041802154-281 / 48 f * 


1 


04 


18 


02 


157 


28 


17 


48 


790 


4.74 


8.69 


-1.10 


4 


75 


V 


Tn^ i onKQ7 i loioni i , 
J u41oUoo ( 4-z!828Ull* 


5 


04 


18 


05 


374 


28 


28 


01 


114 


8.00 


1.67 


0.00 


1 


00 


III 


Tn^ono;i7i i r i7n')o [ ) 
J U4zUz4 / l-f-Z ( [JoooZ 


4 


04 


20 


24 


713 


27 


03 


53 


219 


6. 14 


1.39 


3.94 


1 


00 


III 


J U4z40 ( 4o+z4UzUlD 


4 


04 


24 


57 


485 


24 


02 


01 


687 


0.46 


9.47 


-5.26 


7 


50 


V 


Tn /i Q7 /I n /I 7 i ncininQ 
J U4z f 4U4 f -f-^olUiyo 


. 


04 


27 


40 


473 


26 


10 


19 


340 


4. 18 


3.68 


0.91 


1 


00 


III 


J U4z |4| zl-t-^oUz4UO* 


_ 


04 


27 


47 


210 


26 


02 


40 


531 


1 .23 


10.27 


-4.25 


9 


50 


V 


Tn/ioom 17 i ncioini , 
J U4ZOU14 l+^DlolUl* 


3 


04 


28 


01 


440 


26 


18 


10 


440 


2.00 


1 .02 


2.90 


1 


00 


III 


J 0428 14404- 2b22409 


3 


04 


28 


14 


400 


26 


22 


40 


800 


2.36 


1.41 


2.52 


1 


00 


III 


J U4zoloo44-^t)iiolDy* 


3 


04 


28 


18 


840 


26 


23 


16 


800 


2.41 


1.15 


1.36 


1 


00 


III 


Tn/ioooon i ncnoi 77 , 
JU4zozo4Z4-^DUol / I* 


3 


04 


28 


28 


42 o 


26 


08 


17 


700 


2.65 


3.26 


0.67 


3 


00 


III 


J U4zo44o44-^4z4zUZ 


4 


04 


28 


44 


644 


24 


24 


20 


214 


0.00 


8.89 


-7.19 


6 


25 


V 


J 042*9020 f 4-z!00oUDl* 


3 


04 


29 


02 


074 


25 


58 


06 


132 


0.10 


4.71 


-0.36 


4 


50 


V 


JU4zyizl ( -(-zoloDyz* 


3 


04 


29 


12 


108 


26 


18 


of) 


i 60 


2 . 60 


1.16 


0.82 


1 


00 


TTT 
111 


J U4zy4o44-)-^4oD^z4* 


3 


04 


29 


43 


446 


24 


36 


22 


475 


4.18 


9.01 


0.81 


4 


50 


V 


inion/icci i ni9i /inQ , 
J U4zy4DOl + ^4ol4yo ; f : 


5 


04 


29 


46 


519 


24 


31 


49 


328 


10.48 


1 .48 


0.00 


1 


00 


III 


JU4zy4 f ZO+ZO^y loo 


4 


04 


29 


47 


236 


25 


29 


18 


892 


0.00 


8.43 


-23.70 


8 


00 


V 


Ju4zyODZ^4-^oUoOOy 


4 


04 


29 


56 


228 


25 


03 


55 


904 


2.25 


4.63 


-3.56 


1 


00 


V 


J u4oUUlzU+^04o4oy* 


3 


04 


30 


01 


207 


25 


48 


45 


914 


0.09 


3.60 


1 .49 


3 


75 


V 


Tri/iQni noi i icq nine 

ju4ouiyoi+zooUiyD 


1 


04 


30 


19 


810 


26 


30 


19 


667 


1.56 


6.10 


-5.87 


4 


75 


V 


Ju4oU44yo+ZOOlU/ 1* 


3 


04 


30 


45 


000 


25 


51 


07 


920 


1 .62 


1.63 


1.06 


1 


00 


III 


Tri/iQnoi r: i i/i /inr;o7 , 
JU4olzzl04-Z44yoz I * 


4 


04 


31 




150 


24 


49 


52 


730 


1.56 


5.12 


0.77 


4 


75 


V 


Ju4oliDoU+Z044oz i * 


3 


04 


31 


26 


304 


25 


44 


32 


640 


0.85 


1 .94 


0.78 


1 


00 


III 


J U4olz ( 4U-T"ZOOozol* 


3 


04 


31 


27 


408 


25 


53 


28 


320 


2.22 


0.97 


1.41 


1 


Of) 


III 


J U4oliODO-|-Z41 / lo4* 


5 


n i 

U4 


3 1 


28 


689 


O -1 
Z4 


1 7 
1 / 


18 


44z 


5.34 


5.60 


0.00 


1 


00 


TTT 
ill 


J U4oloozZ-|-Z4iyUlo ; f : 


5 


04 


31 


33 


22!.) 


24 


19 


01 


556 


3.70 


8.23 


0.00 


3 


75 


V 


JU4oloyoD+z4oboo I 


3 


04 


31 


39 


864 


24 


36 


36 


360 


4.44 


1.68 


1.43 


1 


00 


III 


Tn/iQi/tQnn i n^nri/iio , 
J u4ol4oyU+Z4Uy41o ; f : 


5 


04 


31 


43 


901 


24 


09 


41 


857 


5.42 


4.15 


2.09 


1 


00 


III 


Ju4olo4bo4-^4obloo 


3 


04 


31 


54 


672 


24 


36 


16 


200 


7.07 


1.91 


1.97 


1 


00 


III 


JU4olo r 2o+24obloo 


3 


04 


31 


57 


288 


24 


36 


14 


760 


5.97 


2.91 


1 .28 


1 


00 


III 


j u4oioyzi+z4zoiyu* 


5 


04 


31 


59 


219 


24 


23 


19 


021 


2.96 


9.17 


0.00 


3 


75 


V 


JU4ozUzlZ-+-z4lZ4ol* 


■ j 


04 


32 


02 


124 


Z4 


12 


A Q 
4o 


186 


6.46 


4 . 9 8 


0.00 


1 


00 


TTT 
111 


J 0432 1076+2420377* 


5 


04 


32 


10 


764 


24 


20 


37 


788 


6.33 


4.08 


2.94 


1 


Of) 


III 


J U4oz 144 14-^40 bl 1 4* 


4 


04 


32 


14 


4 1 8 


24 


56 


17 


408 


0.00 


6.72 


1.60 


5 


25 


V 


Tfl/IQin C/1 1 1/11 1 on , 

J U4ozzlo4+Z41 Izzl* 


1 


04 


32 


21 


546 


24 


1 1 


22 


106 


0.00 


8.48 


-16.50 


8 


Of) 


V 


JU4ooUiy 1+loolUbb* 


3 


04 


33 


01 


917 


18 


31 


06 


654 


3.01 


3.52 


0.87 


. 


Of) 


III 


TH/1QQ(n77 1 lOT/lnCO, 

JU4ooUz/ f4-l«^4zoo* 


3 


04 


33 


02 


832 


18 


24 


26 


280 


0.34 


1.06 


1.70 


1 


00 


III 


Tn/iQQ/inn7 i noi79')'3 
J04oo420 i+^zl/ ooo 


4 


04 


33 


42 


077 


22 


17 


33 


320 


0.10 


6.84 


1.10 


5 


25 


V 


J U4oo4ol4-f-^D4 / 4z4* 


i 

4 


04 


33 


A Q 
4o 


1 A 1 
14 ( 


26 


1 7 
4 l 


-tz 


■182 


2 . 4 5 


4.22 


-0.52 


4 


2t:> 


V 


TftA o A i /ion i i/iicncn , 
J U4o414zy4-^41bUbz* 


1 


04 


34 


14 


298 


24 


16 


06 


2 1( ) 


2.07 


4.60 


-5.26 


5 


00 


V 


J04o4oo4o-|-22aaU < U* 


_" 


04 


34 


33 


465 




53 


07 


051 


1.17 


8.22 


-4.04 


8 


00 


V 


Tfl^Q/1 CQCC 1 T/imCJO , 

J 04o40ob6+24105oO* 


5 


04 


34 


53 


665 


24 


10 


53 


803 


0.02 


10.40 


-5.73 


6 


25 


V 


J04ooOo2U+22oOo8o 


3 


04 


35 


05 


184 


22 


50 


39 


120 


3.41 


1.39 


0.40 


1 


00 


III 


JU4oul f 044-ZZ4ZZ0Z* 


1 


04 


35 


17 


844 


22 


42 


28 


238 


1.35 


7.61 


-4.99 


4 


75 


V 


J U4oOiy4o+zzoozUo* 


5 


04 


35 


19 


434 


22 


38 


20 


386 


4. 49 


3.83 


0.00 


1 


00 


III 


JU4ooobyu+zb4Uloz 


4 


04 


35 


36 


909 


26 


40 


18 


250 


5.58 


3.19 


0.84 


1 


25 


III 


J 043539 12+22403b9* 


5 


04 


35 


39 


129 


22 


40 


36 


934 


1.13 


6.31 


-7.21 


4 


25 


V 


Tn/iQ^ccnn i i7i Qnnn , 

Ju4ooobyy+^ / louyz* 


4 


04 


35 


57 


000 


27 


13 


09 


250 


0.00 


7. 14 


-3.24 


5 


25 


V 


JU4oDZ ( D0-)-zzoy4yD* 


1 


04 


3 b 


Z I 


653 


22 


39 


4y 


6 5 1 


0.21 


10.67 


-0.72 


6 


■ >f) 


V 


TH/1'371 HOI I i/ionrcjc , 
JU4o i lzzl4-z4oyUoo* 


4 


04 


37 


12 


218 


24 


39 


03 


618 


0.10 


6.84 


-3.06 


5 


75 


V 


J04oolo4U+24ooloO 


4 


04 


38 


18 


405 


24 


55 


13 


505 


0.00 


5.57 


-4.36 


4 


75 


V 


104384651 4-9 5581 53 


3 


04 


38 


46 


512 


25 


58 


14 


880 


3.97 


0.81 


3.72 




00 


III 


J04392692+2552592* 


5 


04 


39 


26 


920 


25 


52 


59 


286 


15.44 


2.70 


0.00 


4 


Of) 


III 


J04393591+2542325* 


2 


04 


39 


35 


916 


25 


42 


32 


598 


0.00 


7.25 


1.12 


6 


75 


V 


J04395172+2549182* 


5 


04 


39 


51 


726 


25 


49 


18 


282 


10.33 


2.82 


0.00 


1 


Of) 


III 


J04401040+2556029 


3 


04 


40 


10 


368 


25 


56 


03 


120 


3.75 


7.12 


-1.47 


4 


75 


V 


J04402392+2554478* 


2 


04 


40 


23 


927 


25 


54 


47 


808 


4.61 


6.74 


-5.05 


5 


25 


V 


J04404085+2542054* 


3 


04 


40 


40 


824 


25 


42 


05 


400 


3.18 


0.99 


3.51 


1 


00 


III 


J04411647+2507537 


3 


04 


41 


16 


464 


25 


07 


53 


760 


2.81 


2.33 


2.58 


1 


00 


III 


J04411897+2515090* 


2 


04 


41 


18 


979 


25 


15 


09 


043 


0.29 


7.63 


-13.08 


9 


25 


V 


J04413108+2519222* 


5 


04 


41 


31 


085 


25 


19 


22 


300 


10.93 


3.26 


0.00 


1 


Of) 


III 


J04422975+2519485 


2 


04 


42 


29 


752 


25 


19 


48 


508 


5.16 


5.95 


-7.44 


4 


75 


V 


J04425194+2414525 


4 


04 


42 


51 


946 


24 


14 


52 


598 


0.51 


4.78 


0.85 


4 


25 


V 


J04430638+2523237* 


1 


04 


43 


06 


382 


25 


23 


23 


766 


3.14 


6.05 


-3.39 


3 


75 


V 


J04430807+2522103* 


3 


04 


43 


08 


040 


25 


22 


10 


560 


3.72 


2.88 


0.38 


1 


Of) 


III 


J04435916+2419389* 


4 


04 


43 


59 


169 


24 


19 


38 


989 


1.09 


4.62 


-3.81 


4 


75 


V 










Uncertain mcmbci 


ship 












.704272297+2636474 


4 


04 


27 


22 


975 


26 


36 


47 


498 


2.41 


5.49 


-6.94 


5 


50 


IV 


J04315129+2506524* 


4 


04 


31 


51 


298 


25 


06 


52 


488 


1.86 


3.94 


-2.62 


4 


75 


IV 


J04355760+2253574* 


5 


04 


35 


57 


608 


22 


53 


57 


491 


2.55 


5.57 


-21.58 


4 


75 


IV 



Table 3. Derived parameters and spectral properties of the sources identified as field stars or with uncertain Taurus 
membership. The first column gives the 2MASS name. The second column corresponds to the observation run as indi- 
cated in Table |2 The following two columns are 2MASS coordinates. Reddening values, spectral type and luminosity 
class derived from the spectral fitting procedure (see t!3.1.1[) are listed in columns (5), (8) and (9) respectively. Columns 
(6) and (7) give the measured Na I and Ha equivalent widths. The 2MASS names followed with a * correspond to 
candidates not selected with the primary selection method (see H2.1|) . 




Fig. 5. Medium-resolution optical spectra of the new Taurus very low-mass members and brown dwarfs. The spectra 
have been corrected for extinction. The derived spectral type is indicated. The spectrum of ZZ Tau IRS is also shown. 
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2MASS Name 


Name 


Run 


RA(J2000) 


DEC(,J2000) 


Av 


EW(Na) 


EW(H„) 


SpT 




J04151471+2800096 


KPNO-Tau-1 


1 


04:15:14.714 


28:00:09.612 


0.41 


3.40 


-26.82 


9.00 


IV 


J04174955+2813318 


KPNO-Tau 10 


3 


04:17:49.554 


28:13:31.854 


0.00 


2.78 


-40.09 


5.50 


IV 


J04185115+2814332 


KPNO-Tau-2 


5 


04:18:51.156 


28:14:33.241 


0.37 


3.16 


-8.39 


6.75 


IV 


J04210795+2702204 


CFHT-Tau 19 


4 


04:21:07.953 


27:02:20.418 


7.34 


2.17 


-442.17 


5.25 


IV 


J04214631+2659296 


CFHT-Tau 10 


4 


04:21:46.315 


26:59:29.609 


3.59 


3.49 


-17.32 


6.25 


IV 


J04221644+2549118 


CFHT-Tau 14 


4 


04:22:16.440 


25:49:11.842 


0.57 


3.34 


-48.77 


7.75 


IV 


J04221675+2654570 


CFHT-Tau 21 


4 


04:22:16.759 


26:54:57.078 


6.61 


2.25 


-45.09 


1.25 


IV 


J04242646+2649503 


CFHT-Tau 9 


4 


04:24:26.462 


26:49:50.362 


0.91 


3.39 


-9.95 


6.25 


IV 


J04272799+2612052 


KPNO-Tau-4 


2 


04:27:27.997 


26:12:05.270 


2.45 


2.07 


-158.08 


9.50 


IV 


J04274538+2357243 


CFHT-Tau 15 


4 


04:27:45.380 


23:57:24.325 


1.30 


3.76 


-18.90 


8.25 


IV 


J04292165+2701259 


CFHT-Tau 18 


4 


04:29:21.653 


27:01:25.946 


4.85 


1.77 


-38.30 


6.00 


IV 


J04295950+2433078 


CFHT-Tau 20 


4 


04:29:59.508 


24:33:07.852 


3.60 


2.00 


-114.03 


5.50 


IV 


J04300724+2608207 


KPNO-Tau-6 


2 


04:30:07.244 


26:08:20.792 


0.88 


3.51 


-207.91 


9.00 


IV 


J04302365+2359129 


CFHT-Tau 16 


4 


04:30:23.655 


23:59:12.991 


1.51 


4.05 


-16.76 


8.50 


IV 


J04305171+2441475 


ZZ Tau IRS 


4 


04:30:51.714 


24:41:47.512 


2.38 


2.25 


-141.13 


5.25 


IV 


J04312669+2703188 


CFHT-Tau 13 


4 


04:31:26.690 


27:03:18.810 


3.49 


4.20 


-3.60 


7.25 


IV 


J04321786+2422149 


CFHT-Tau 7 


3 


04:32:17.862 


24:22:14.984 


0.00 


2.83 


-8.63 


6.50 


IV 


J04325026+2422115 


CFHT-Tau 5 


2 


04:32:50.265 


24:22:11.564 


9.22 


3.98 


-29.84 


7.50 


IV 


J04330945+2246487 


CFHT-Tau 12 


4 


04:33:09.457 


22:46:48.702 


3.44 


2.51 


-79.65 


6.50 


IV 


J04350850+2311398 


CFHT-Tau 11 


4 


04:35:08.508 


23:11:39.865 


0.00 


3.13 


-45.07 


6.75 


IV 


J04390396+2544264 


CFHT-Tau 6 


3 


04:39:03.960 


25:44:26.416 


0.41 


3.29 


-63.74 


7.25 


IV 


J04400174+2556292 


CFHT-Tau 17 


1 


04:40:01.745 


25:56:29.227 


6.50 


2.87 


-7.26 


5.75 


IV 


.104411078+2555116 


CFHT-Tau 8 


3 


04:41:10.785 


25:55:11.651 


1.77 


2.89 


-52.00 


6.50 


IV 



Table 4. Same as Table for confirmed Taurus members. 



of 1.17 (see Figure^}. However a foreground nature does 
not seem compatible with the near-infrared colors. It is 
therefore difficult at this stage to conclude on the nature 
of this source. 

The remaining two sources, 

2MASS J04272297+2636474 and 

2MASS J04315129+2506524, have mid-M spectral 
types, low-level line emission (EW(Hq) emission <8 A) 
and occupy in Figure 01 the same location as dwarfs. For 
these two sources, spectral types M5.5 and M4.75 and 
Ay values of 2.41 and 1.86 are derived from the best fit 
solution obtained with the average of a dwarf and giant 
template. Their derived location in the HR diagram with 
the above parameters is shown in Figure [3 Both sources 
fall between the 10 and 30 Myr isochrones, at the low 
luminosity end of the Taurus population. However, as 
illustrated in Figure the fits obtained with a dwarf 
template alone appear in both cases very close to the 
previous best fit solution. The differences in x 2 between 
the dwarf and sub-giant fits for these two sources are 
0.06 and 0.35 respectively, significantly lower than the 
average of the x 2 differences for the remaining Taurus 
candidate population of 1.8±1.25. The dwarf fits give a 
similar spectral type but significantly lower reddening 
values compatible, within our derived 1 a uncertainty of 
0.8 mag, with a foreground field dwarf nature. Their near- 
infrared colors ((J-H)=0.35, (H-K) =0.7-0.8) and I band 
magnitudes appear also compatible with a foreground 
nature. 

We therefore classify these last 3 sources in Table |3 
with uncertain membership and do not consider them as 
Taurus members in the following. We arc therefore left 
with 23 bonafide Taurus members. We note that all but 2 
of these show strong Ha emission with equivalent widths 
in excess of 8 A. All these 23 sources were included in our 
primary sample of 51 sources selected according to the 
criteria of § 12.11 We therefore estimate the efficiency of 
our selection procedure to be 45 %. 



4.2. Properties of the new Taurus members 

We summarize in Table |S] the derived properties of the 
Taurus members. We list the spectral types and red- 
dening values derived through our spectral fitting pro- 
cedure. The effective temperature is estimated from the 
spectral type using the following linear fit to the rela- 
tionship derived bv iLuhmanl (^.999) for young M dwarfs: 
T eff = 3780 - 175.6 * SpT. We compute the luminosity 
from the dereddened I band magnit ude, using the bolo- 
metric corrections for M dwarfs from lLuhmanl 1 1999ft . We 
also list our / and z' band photometry combined with the 
J, H and K 2MASS photometry. We show in Figure|7|thc 
derived HR diagram for our new Taurus sources, includ- 
ing previously known members . We u se the 0.08 M Q mass 
track from the IChabrier et alJ l)2000h pre-main sequence 
models to define the stellar/substellar boundary, which 
falls at spectral type M6.25V for a typical Taurus age of 
3 Myr. We will therefore use in the following this spectral 
type as the limiting boundary for brown dwarfs in our 
sample. We note in the substellar domain a general trend 
of the models to predict larger ages for the lowest mass 
objects. This tendency is likely due to uncertainties in the 
pre-main sequence tracks and/or bolometric corrections 
for the lowest masses. 

We recover in our sample 6 previously known Taurus 
members, and we find 12 new BDs and 5 new very low 
mass (VLM) stars. We show in Figure [3] the individ- 
ual spectra of the 17 new Taurus members. We also 
show ZZ Tau IRS, for which, to our knowledge, no low- 
resolution optical spectrum is published in the literature. 
Four of the new BD s were previously pu blish ed in confer- 
ence p roceedings bv lMonin et alJ l)2004|) and lGuieu et alJ 
l(2004h . Two of these, CFHT-Tau 6 (J04390396+2544264) 
and CFHT-Tau 8 (J0441 1 078+25551 16), have indepen- 
dently been discovered bv ILuhmanl l)2004h . This author 
found the same spectral type for CFHT-Tau 6 (M7.25) and 
a close extinction value (Ay=0.26). For CFHT-Tau 8 we 
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X (A) X (A) X (A) 

Fig. 6. Spectral fits (red thick lines) obtained for the 3 sources (black thin lines) with uncertain membership. The 
top panels represent the best solution obtained with a dwarf template while the bottom panels show the best fits 
obtained with the average of a dwarf and giant. The vertical dashed lines represent the wavelength limits used for the 
fit (7000-8500 A). 



Name 


ST 


Av 


logi (L/Lsun) 


Teff 


I 


I-z 


J-H 


H-K 


K 


ZZ Tau IRS 


M5.25 


2.4 


-1.4627 


3100 


15.95 


0.85 


1.40 


1.13 


10.31 


KPNO-Tau 10 


M5.50 


0.0 


-1.3214 


3065 


14.11 


0.11 


0.77 


0.32 


10.78 


KPNO-Tau-2 


M6.75 


0.4 


-2.1302 


2889 


16.56 


0.49 


0.68 


0.49 


12.75 


KPNO-Tau- 1 


M9.00 


0.4 


-2.6340 


2571 


18.15 


0.54 


0.85 


0.47 


13.77 


KPNO-Tau-6 


M9.00 


0.9 


-2.5129 


2571 


17.90 


0.66 


0.80 


0.51 


13.69 


KPNO-Tau-4 


M9.50 


2.5 


-2.3744 


2500 


18.75 


0.71 


0.97 


0.74 


13.28 



CFHT- 


Tau 


21 


Ml. 25 


6.6 


-0.4185 


3665 


15.17 


0.80 


1.54 


1.03 


9.01 


CFHT- 


Tau 


19 


M5.25 


7.3 


-1.1423 


3100 


16.54 


1.15 


1.79 


1.52 


10.54 


CFHT- 


Tau 


20 


M5.50 


3.6 


-0.8586 


3065 


15.10 


0.79 


1.14 


0.73 


9.81 


CFHT- 


Tau 


17 


M5.75 


6.5 


-1.1652 


3030 


17.62 


0.82 


1.58 


0.88 


10.76 


CFHT- 


Tau 


18 


M6.00 


4.9 


-0.3715 


2994 


14.63 


0.97 


1.30 


0.77 


8.73 


CFHT- 


Tau 


10 


M6.25 


3.6 


-1.6792 


2959 


16.85 


0.99 


1.09 


0.60 


12.13 


CFHT- 


Tau 


9 


M6.25 


0.9 


-1.6218 


2959 


15.35 


0.77 


0.69 


0.43 


11.76 


CFHT- 


Tau 


7 


M6.50 


0.0 


-1.2149 


2924 


14.12 




0.73 


0.39 


10.40 


CFHT- 


Tau 


8 


M6.50 


1.8 


-1.6940 


2924 


16.43 




1.05 


0.68 


11.45 


CFHT- 


Tau 


12 


M6.50 


3.4 


-1.4401 


2924 


16.26 


1.03 


1.01 


0.59 


11.55 


CFHT- 


Tau 


11 


M6.75 


0.0 


-1.5964 


2889 


14.88 


0.71 


0.59 


0.35 


11.59 


CFHT- 


Tau 


13 


M7.25 


3.5 


-2.1012 


2818 


17.90 


1.05 


0.86 


0.52 


13.45 


CFHT- 


Tau 


6 


M7.25 


0.4 


-1.6223 


2818 


15.40 


0.53 


0.80 


0.47 


11.37 


CFHT- 


Tau 


5 


M7.50 


9.2 


-1.1258 


2783 


18.79 


1.09 


1.74 


0.94 


11.28 


CFHT- 


Tau 


14 


M7.75 


0.6 


-1.7494 


2747 


15.62 


0.88 


0.70 


0.42 


11.94 


CFHT- 


Tau 


15 


M8.25 


1.3 


-2.4380 


2677 


17.94 


1.08 


0.69 


0.55 


13.69 


CFHT- 


Tau 


16 


M8.50 


1.5 


-2.4152 


2642 


17.91 


1.07 


0.72 


0.54 


13.70 



Table 5. Photometry and derived parameters for previously known (top) and new (bottom) Taurus BD and VLM 
members 



derive an M6.5 spectral type wi th a moderate reddening 
(A v = 1.8) while iLuhmanl l)2004h found an M5.5 spectral 
type with Aj=0.7 (Av=2.6). We note that the discrep- 
ancy in Ay values is similar to our estimated reddening 
uncertainty (0.8 mag) while the difference of 1 spectral 
class is twice our estimated spectral typing error. 



Among the remaining Taurus members, our estimates 
of spectral type and reddening values appear generally 
compatible, within the quoted uncertainties, with values 
given in the literature, except for 2 sources for which larger 
discrepancies are found: ZZ Tau IRS and KPNO Tau 4. 
ZZ Tau IRS, the suspected energy source of HH 393 
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ijGomez et al.lll997l). has been recently ob served spectro- 
scopicaliv bv I White & Hillenbrand! <)2004[K These authors 
derive from their high-resolution optical spectroscopic ob- 
servations an M4.5±2 spectral type, compatible within the 
quoted uncertainties with our own derivation of M5.25. 
However, using the (J-H) color, they estimate a reddening 
value of 7.6, significantly larger than our own estimate of 
Ay = 2.4, derived from the spectral fitting procedure. We 
note that this source displays a significant excess in the ( J- 
H)/(H-K) diagram as well as strong accretion signatures, 
which may significantly affect the derivation of redden- 
ing from near-infrared photometry. Indeed, if we deredden 
ZZ Tau IRS back to the locus of accreting T Tauri stars 
in the (J — H)/(H — K) diagram, we find Ay = 3.1, a 
value compatible with our derivation from fitting of the 
optical low resolution spectrum. We note however that 
these values may be a lower limit to the true redden- 
ing if indeed this source is viewed edge-on as suggested 
by its unusually l arge emission lines equ i valent widths, as 
already noted bv IWhite fc Hillenbrand! l)2004|) . We come 
back to this issue in the following section. In the case 
of KPNO Tau 4, we find the same spectral type (M9.5) 
but a reddening val ue larger than the one computed by 
iBriceno et all (|2002h (A v = 2.5 vs 0.0). This discrepancy 
likely arises from the rapid increase of the continuum flux 
longwards of ~ 8200 A for spectral types later than M9V. 
This part of the continuum therefore strongly influences 
the derivation of reddening values at late spectral types. 
Uncertainties on its exact shape result in large uncertain- 
ties on derived Ay. We therefore estimate that for our 
later spectral types (> A/9), uncertainties on derived Ay 
values are almost doubled and likely increase to 1.6 mag. 




3.6 3.5 3.4 



log(Teff) 

Fig. 7. Derived H-R diagram in Taurus. New Taurus 
members identified in this work are indicated with Star 
symbols. Previously known Taurus members are shown 
as open circles. We have also plotted with triangles the 
3 sources with uncertain members hip discussed in § 4.1 . 
The pre-main sequence tracks from lChabrier et alJ l|2000F ) 
(DUSTY models) are shown as blue curves: solid lines rep- 
resent isochrones with ages of 1, 3, 10, 30 and 100 Myr, 
dotted lines represent mass tracks ranging from 0.02 M Q 
to 1.4 M Q as labelled. The green curve shows the com- 
puted stellar /substellar boundary at 0.08 M Q . 



4.3. Accretion and outflow signatures 

We examine the Taurus VLM members identified in this 
work for evidence of accretion/outflow. These signatures 
include strong near-infrared excess, suggestive of an ac- 
tive accretion disk; strong emission in the optical forbid- 
den lines which, in Classical T Tauri stars, trace out- 
flow signatures strongly correlated with the a ccretion pro- 
cess l|Hartigan et al.lll994tlCabrit et al.lll990|) : strong Ha 
emission beyond the limit for chromosphcric activity. A 
large fraction of the newly identified Taurus VLM mem- 
bers show indications of one or more of these accre- 
tion/outflow signatures. 

Three of the VLM Taurus members (CFHT-Tau 19, 
ZZ Tau IRS, CFHT-Tau 21 ) and one of the BDs (KPNO- 
Tau 4) exhibit near-infrared excess in the (J-H) versus 
(H-K) diagram (Figure |HJ. These sources lie in this di- 
agram to the right of the reddening vector extending 
from the M7 dwarf color and above the locus of accret- 
ing T Tauri stars, suggesting near-infrared excess emis- 
sion arising from the presence of an accretion disk. Using 
our estimates of reddening and spectral type, we further 
compute a A(H-K) excess for each confirmed Taurus mem- 
ber. This near- infrared excess is represented in Figure ITUI 



as a function of spectral type and compared to the maxi- 
mum (H-K) color excesses predicted for a flat reprocessing 
circumstellar disk around a c entral star of sp ectral type 
ranging from M0 V to M9 V l|Liu et al.ll2003j) . The three 
VLM members previously identified in the (J-H) / (H-K) 
diagram show a clear (H-K) color excess above what is 
expected from a pure reprocessing disk. Given the one a 
uncertainty of 0.1 mag in A (H-K), we do not find evi- 
dence from this diagram for strong active accretion in the 
substellar population. However, the lack of detectable (H- 
K) excess does not imply an absence of disk around these 
sources. Indeed, detailed modelling of the spectral energy 
distributions from reprocessing disks a round young sub- 
stellar objects by iNatta fc Test] l(200ll) showed that ex- 
cess disk emission becomes clearly detectable longward of 
3 /im only. 

The three sources with strong near-IR excess also dis- 
play strong emission in the optical forbidden lines of 
[O i] A6300, 6363 A, [N n] A6548, 6583 A, and [S n] A6716, 
6731 A (Figure EJ, which, in higher mass CTTs, trace 
outflow signa tures strongly correlated with the accre- 
tion process l(Hartigan et all Il994t ICabrit et"afl Il990h . 
One of these sources, ZZ Tau IRS, shows anomalously 
large equivalent widths in the optical forbidden lines 




S. Guieu et al.: 17 new very low-mass members in Taurus 

40 



13 



(H-K) 

Fig. 8. ( J-H) / (H-K) diagram. VLM stars (open star sym- 
bols) and brown dwarfs (filled star symbols) Taurus mem- 
bers identified in this work are represented. Numbers iden- 
tify the 4 VLM stars with intense line emission from 
FigureEl Also shown are the locations of the giant branch 
(open green triangles) and the dwarf sequence (thick black 
curve) extending from G2 to M7 (SpT). The black thin 
line illustrates a reddening vector of 10 Ay extending from 
the M7 V dwarf co lors. The locus of th e accreting classical 
T Tauri stars from lMever et all 1 1997ft is displayed with a 
dotted-dashed line. 

(EW([0 i] A6300 A) = -190 A, EW([S II] A6731+6716 A) 
= -77 A). Typical equivalent widths in T Tauri stars are 
on the order of a few tenths of A. Such large values in 
ZZ Tau IRS could be explained if the source is seen close 
to edge-on and the continuum is therefore suppressed rel- 
ative to the line emission, arising from extended jet emis- 
sion. The suggestion that Z Z Tau IRS is close to edge-on 
has already been made by IWhite k, Hillenbrand! l|2004|) 
on the same argument. One VLM object , CFHT-Tau 20 
(J04295950+2433078) has signific ant Ca II line emission , 
an indi cator of accretion acti vity l|Muzerolle et alJll998l) . 
Indeed. iMohantv et alJ l)200fj) demonstrated recently that 
the Can A8662 A line flux correlates very well with the 
disk mass accretion rate down to the substellar regime. 

As noted before, all but 2 of the Taurus members iden- 
tified in this study display strong Ha emission with equiv- 
alent widths in excess of 8 A. This level of Ha emission 
is usually taken in higher mass T Tauri stars (TTs) as 
a strong accretion indicator. Chrosmospheric Ha emis- 
sion activity is however expected to increase strongly 
towar ds late M types. iBarrado v Navascues fc Martini 
l)2003(l conducted an extensive study of Ha emission in 
a sample of young M stars. They derived an empirical 
CTTs/WTTs boundary, extended to substellar analogs, 
in the Ha equivalent width versus spectral type dia- 
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Fig. 9. A zoom on representative regions of the spectra is 
presented for the 4 Taurus members with strong emission 
lines. 
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Fig. 10. (H-K) color excess versus M spectral type for 
Taurus members from this study. VLM stars with intense 
line emission are numbered as in Figure |U| The solid line 
shows the maximum (H-K) color excess expected from a 
flat reprocessin g circumstellar disk around a MOV to M9V 
central source iLiu et al.ll2003T) . 
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Fig. 11. Ha equivalent width versus spectral type for 
VLM and BD Taurus members. Star symbols identify 
the Taurus members studied in this paper. Previously 
known members are shown with black dots. The dashed 
line shows the empirical boundary between chromo- 
spheric and accretion-related Ha a ctivity derived by 
iBarrado v Navascues fe Martini l|2003|) . The EW {Ha) is 
plotted on a logarithmic scale. 



gram. This boundary is defined by the saturation limit 
for chromospheric activity at log (L# Q /Lboi) = —3.3, de- 
termined from study of a sample of young open clus- 
ters. We plot in Figure ^2 the distribution of late-type 
Taurus members with respect to this boundary. We in- 
clude members from this study and previously known 
late-type sources with published Ha equivalent width s 
measurements fromlLuhmanl ll2000l);lMartrn et akl l|200l|) ; 
IBricefio et all ll2002h : lLuhma,n et al.l (|2003ah. Figure El 
suggests that 6 Taurus brown dwarfs among our sample 
( CFHT-Tau 6, CFHT-Tau 11, CFHT-Tau 12, CFHT- 
Tau 8, KPNO-Tau-4 and KPNO-Tau-6) have level of 
Ha emission in excess of chromospheric activity and 
are likely CTTs substellar analogs. All but one (CFHT- 
Tau 17) of the VLM members from this study lie above 
the CTTS/WTTs boundary. In particular, the 4 VLM 
sources with strong emission lines (labelled 1 to 4), dis- 
cussed in the previous paragraph, clearly stand out in this 
diagram as strongly accreting sources. We note that high- 
resolution spectroscopic observations that fully resolve the 
Ha profiles are required to confirm that accretion is in- 
deed present in these sources. 

The near-infrared excess emission and line emission 
properties of the Taurus VLM members analyzed here 
are strikingly similar to the accretion/ejection signatures 
observed in the higher mass CTTs, indicating that the 
same processes are occurring down to the substellar limit. 
The apparent large fraction of accreting sources in our 
VLM sample is probably a result of our selection proce- 



dure since our sample of objects with spectral types earlier 
than M4 V is biased towards large near-IR excess sources. 
The fraction of accreting brown dwarfs, as deduced from 
the properties of their Ha emission, is approximatively 
40% (6/15) in our sample and 42 % (10/24) when we in- 
clude previously known substellar members. We note that 
the substellar sample studied here does not suffer from 
the same bias as the higher mass objects and represents 
72 % of all currently known BDs in Taurus. This derived 
substellar accretor fraction in Taurus is statist i cally con- 
sistent with the one derived bv lMohantv et alJ l|2005h for 
sources with spectral type > M5 from a detailed analysis 
of accretion signatures of 32 % ± 10 %. The fraction of 
accreting sources among the substellar population is simi- 
lar to the one in the higher mass T Tauri stars, suggesting 
that inner disk lifetimes around young Taurus BDs appear 
on the same order or larger than around higher mass TTs. 



4.4. The stellar to brown dwarf ratio revisited 

Following IBricefio et all (|2002jk iLuhman et all l)2003a|) : 
Luhman ( 2004) , we quantify the relative numbers of brown 
dwarfs and stars in Taurus. The completeness in mass 
of our spectroscopic survey is set on one hand by the 
completeness limits of the 2MASS observations (J=15.25, 
H=14.4, K=13.75), which have been combined with the 
optical data to select substellar Taurus candidates, on the 
other hand by the i' band limit of 20 of our spectroscopic 
observations. Using t he pre-main sequence tracks from 
IChabrier eta l. (2000) (DUSTY models), we compute in 
both cases a corresponding mass limit of 30 Mj and red- 
dening limit of Ay < 4 for a maximum age of lOMyr at 
the Taurus distance (140 pc). 

We derive below the reddening and mass limited sam- 
ple of all known Taurus members and update the current 
estimate of the brown dwarf to star ratio by including the 
new fields studied in this paper. 

Co mpiling the list of kn own Tau rus members 
from IKenvon fc Hartmannl <ll995h : | Bricer1o et al 



(1998, 
(2000) 



iKenyon fc Hartmannl lll995lk | Brice 
I1999L 120021): ILuhman et all lll998h; [ 




.Luhmarj 

Luhman et al J l | 2003al): ILuhman! 

[ White fc Hillenbrand! ll2004h: IStrom fc Stroml , 

iMartm et alJ l|200ft l200ll) : IWichmann et all l|2000l) and 

including the 17 new members from this work, we have 
identified 64 Taurus members (including known compan- 
ions) projected towards o ur new fields but no t included 
in the c ombined surveys of IBricefio et alJ Jl998h : lLuhmar] 
d2000l): IBricefio et alJ l|2002l) : ILuhm an! l(2004h Following 
IBricefio et al.1 (|2002h . binary systems with separations 
< 2" are considered as individual objects. We are therefore 
left with 54 sources. IWhite fc Hillenbrand! l)2004l) recently 
provided estimates of spectral types and reddenings 
from high-resolution optical spectroscopy for a sample of 
embedded young stars in Taurus. We adopt their derived 
effective temperatures, luminosities and Ay values for 
Haro 6-5B and 5 embedded IRAS sources included in 
our sample (IRAS 04016+2610, IRAS 04260+2642, 
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IRAS 04264+2433, IRAS 04278+2253A, 
IRAS 04303+2240). We are left with 11 embedded IRAS 
sources for which we currently la ck accurate estimates 
of spe ctral types and reddenings. IWhite fe Hillenbrand! 
l)2004() derive, for the embedded Class I sources in Taurus 
for which they have accurate spectral types, an average 
visual extinction of 10.5 magnitudes with a standard 
deviation of 4.2 magnitudes. The 11 IRAS sources are 
therefore also likely to be very embedded and not to 
contribute significantly in numbers to our reddening lim- 
ited sample. As a consistency check, we derive for these 
sources a rough estimate of Ay from the (J-H)/(H-K) 
diagram following the method described in i )3.1.2l Among 
the 11 IRAS sources, only one source ends up with 
Ay < 4. We further exclude F T Tau, a strong continuum 
source ilHerbig fc Belli Il988l). and IRAS 04016+2610 
and ZZ Tau IRS for which IWhite fc Hillenbrand! l|2004 
present strong indications that these are edge-on systems. 
Among the remaining 40 sources, 31 sources have Ay < 4 
and are therefore included in the reddening limited 
sample, 9 of these arc BDs. 

J04202555+2700355, J04213459+2701388, 
J04284263+27 14039 a n d J0 442713+2512164, sources 
discovered by iLuhmanl l)2004h . are not included in the 
substellar/stellar ratio computed by this author because 
they correspond to an initial and incomplete 225 deg 2 
survey. However, since they are also present in our 
fields, they are included into our present estimate of the 
substellar to stellar ratio. In addition, the new Taurus 
members, CFHT-Tau 17, CFHT-Tau 5, CFHT-Tau 7 
and CF HT-Tau 1 1 disco v ered in this study, fall i nside the 
fields of ILuhmanl l|2000|) : IBriceno et all (|2002j) : ILuhmanl 
(200^) but are not reported by these authors. CFHT- 
Tau 17 and CFHT-Tau 5 are embedded (Ay = 6.5 and 
Ay = 9.2) but CFHT-Tau 7 and CFHT-Tau 11 have 
little reddening and are therefore included in our ratio. 
We estimate a spectra l type of M6.5 for CFHT-Tau 8 
different from iLuhmal l(2004h (M5.5), so this source falls 
in our sample under the hydrogen burning limit. We note 
that CFHT-Tau 7 lies at a projected distance of only 20" 
from V928 Tau A/B and could be a wide companion to 
this binary system. 

We combine the additional reddening limite d sample 
derive d from this study with the one defined bv ILuhmanl 
l)2004h and compute an updated substellar to stellar ratio 
in Taurus. We first use the same mass completeness limit 
of 20 Mj as previous studies for comparative purposes: 
K ss = 7V(0.08 - O.O2Af )/7V(O.O8 - 1OM ) = 31/127 = 
0.24 ± 0.05. Using a mass completeness limit of 30M,/ in- 
stead results in K ss = 29/12 7 = 0.23 ± 0.05 

As already mentioned in ILuhmanl l|2004h . the derived 
TZ SS ratio in Taurus is likely an upper limit for possible 
remaining incompleteness at low stellar masses. Indeed, 
the distribution of spectral types in Tauru s shows an ap- 
parent deficit at spectral types M2 to M4 l|Luhman et alJ 
20033). Searches for substellar objects in Taurus are, as 
this study, typically complete for spectral types later than 
M4, while previous optical and X-ray large-scale surveys of 
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Fig. 13. Normalized histograms of the distance to the 
nearest stellar neighbor, for stellar members only (filled 
histogram) and for the substellar population only (empty 
histogram) . 



the Taurus cloud are likely complete down to ~ M2. The 
increase of Rss with increasing spatial coverage in Taurus 
could be due to a larger incompleteness at these early to 
mid-M spectral types in the distributed population than in 
the more studied high stellar density regions. However, we 
show in the next section that this is likely not the case: the 
distributed population does not show an apparent larger 
deficit of M2-M4 stars (with respect to the earlier spectral 
type population). Moreover, the incompleteness at spec- 
tral types M2-M4 in Taurus has a limited effect on the 
computed Rss ratio: increasing by a factor 2 the known 
number of stars in this spectral range would bring the Rss 
ratio down to 0.20 +/- 0.05, compatible within one sigma 
with our new estimate. 

Our new estimate of the Taurus substellar to stellar 
ratio of 0.23 +/- 0.05 is higher than the value derived by 
Luhman (2004) of 0.18 +/- 0.04, although still compat- 
ible at the 1 sigma level. It is now however comparable 
to the Trape zium value of 7^ s , s (Tr apezium) = 0.26 ± 0.04 
estimated by IBriceno et alJ l)2002l) , using the same evolu- 
tionary models and treating binary systems in the same 
manner. Recent estimates seem to sugge s t an e ven lower 
value for the Trapezium. ISlesnick et alJ l|2004h have re- 
cently computed a substellar mass function for the Orion 
nebula cluster with a sample of spectroscopically con- 
firmed members. These authors, using a differen t pre-main 
sequence model ( D'Antona fc Mazzitellil ll997tL found a 
ratio of 7£ ss =0.20, lower than the pre vious estimates 
derived from photometric studies alone jLuhman [2000; 
Iffillenbrand fc Carpenterll200d: iMuench et alJl2002|) . 
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Fig. 14. Variation of the ratio of substellar to stellar mem- 
bers in the aggregates (black dots) versus the distributed 
population (red triangles) as a function of the radius taken 
for the aggregates (i n units of the physical radius R c from 
iGomez et al.l l|l993fl 'l. The vertical arrow corresponds to 
an equal number of members in the aggregates and dis- 
tributed population. The curve below shows the corre- 
sponding average stellar density per pc~ 2 within the ag- 
gregates. 



4.5. The spatial distribution of BDs 

Figure El shows the spatial distribution of the cur- 
rently identified stellar and substellar members in Taurus, 
including all the new sources identified in this work. 
Figure El suggests that the Taurus substellar members 
closely follow the distribution of stars. In order to quan- 
tify this effect, we compute the distribution of the near- 
est stellar neighbor, both for the stellar sample (defined 
as members with spectral types earlier than M6.25V) and 
the substellar sample. In this calculation, we consider only 
sources projected towards fields searched for substellar 
members and all binary systems with separation less than 
10" (ie including virtually all known binary systems in 
Taurus) as a single stellar source. This results in 184 stel- 
lar systems and 33 substellar members. The histogram of 
nearest stellar neighbor distances is shown in Figure IHfl 
for both populations. The two distributions appear statis- 
tically consistent. The median of the distribution for the 
stellar sample is: 4.82' with a standard deviation of 8.45', 
while for the substellar sample we find a median of 5.72' 
(standard deviation 10.48'). Considering only the sources 
with Ay < 4 does not significantly change the result. 

However, the fact that the substellar to stellar ratio 
(1Z SS ) has been steadily increasing as the area covered in 
Taurus became larger suggests a dependency of this ratio 
on the stellar density. To quantify this effect, we divide the 
Taurus population in two components: aggregates and dis- 
tributed population. To define the former population, we 



and a fixed factor of the physical radii R c determined for 
each aggregate by these authors, which range from 0.5 to 
1.1 pc. We plot in FigurelTHthe variation of the substellar 
to stellar ratio in both populations as the surface com- 
prising the aggregates is increased. In the computation 
of the 1Z SS ratio we consider only sources with Ay < 4 
and binary systems with separations > 2" are treated as 
a single source. Figure El clearly shows a marked increase 
of the 1Z SS ratio as the integration radii defining the ag- 
gregates are increased. For integration radii ~ 0.7R C , the 
Tt S8 ratio is a factor 2.3 lower in the aggregates than in 
the remaining distributed population. At integration radii 
> R c both 1Z SS values converge towards the average value 
of 0.24 determined earl ier. Figu re l*HH ill ustr ates how 1Z SS 
increases from ~ 0.13 IjBricefio et al.l 12002) towards the 
higher stellar density aggregates, to 0.24 (this work) as 
larger and lower stellar surface density areas are being 
surveyed. Figure El suggests a BD deficit (of a factor ~ 
2) in the central regions of the Taurus aggregates, on spa- 
tial scales ~ 0.35 to 0.7 pc, where stellar densities are the 
highest. 

As mentioned above, this apparent increase of 1Z SS 
could be due to a larger incompleteness in spectral types 
M2V to M4V in the distributed population than in the 
more studied aggregates. To test this hypothesis, we have 
computed the evolution of the ratio N ^gpx<Ai2)^ • ^ e 
find that this ratio is lower in the cores of the aggregates 
than in the distributed population and increases as the 
average stellar density decreases. We conclude that there 
is no strong reason to believe that the incompleteness in- 
creases away from the aggregates and that the observed 
increase of the 1Z SS ratio with decreasing average stellar 
density in Taurus is real. 

5. Discussion: implications for substellar formation 

Two main classes of models have been proposed for the 
formation of substellar objects. In the standard formation 
scenario brown dwarfs form like stars, through (turbu- 
lent) gravitational collapse and fragmentation of cores, fol- 
lowed by subsequent disk accretion. In the ejection models, 
brown dwarfs are stellar embryos ejected from their parent 
core either early in their evolution fr om dynamically un- 
stabl e multiple protostellar systems 1 Reinurth fc Clarkel 
l200l|) or thro ugh secular dynamical decay in dense embed- 
ded cl usters l|Sterzik fc Durisenll2003t iKroupa fc Bouvierl 
l2003al) . Thus deprived of a significant mass reservoir, these 
ejected embryos end up as very low mass objects. We dis- 
cuss below the implications of our observational results 
on the substellar population of the Taurus cloud for these 
two main classes of models. 

The fact that the abundance of BDs (down to 30 Mj) 
relative to stars is found to be the same (~ 25 %) in the dif- 
fuse Taurus cloud and the high-density ONC cluster seems 
to suggest that there is no strong dependency of the sub- 
stellar IMF on initial molecular cloud conditions, in partic- 
ular gas density and level of turbulence. This result seems 
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to be in contrast with the most rece nt simulations of su- 
person ic turbulent fragmentation bvlPadoan fe Nordlundl 
i|2004fk iDekado-Donate et all l|2004|) . with sonic rms 
Mach numbers in the range 3 to 20, which predict a strong 
dependency of the abundance of BDs with the initial cloud 
velocity dispers i on and gas density. In the subsonic case, 
lOoodwin et alJ l(2004bl) show that 20 % of the objects 
formed are brown dwarfs almost independently of the level 
of turbulence. It is unlikely however that a subsonic model 
would apply to the ONC precursor where supersonic tur- 
bulent motions were likely present. It should be noted that 
turbulent fragmentation models predict a decrease of the 
average mass formed with an increasing level of turbu- 
lence, which could account for the difference in peaks in 
the high mass end of the mass distributions of the Taurus 
and ONC populations. 

Although we do not find statistical evidence for a spa- 
tial segregation between stars and brown dwarfs in Taurus, 
we do find strong indication that the abundance of BDs 
relative to stars in the centers of the Taurus aggregates, 
on scales ~ 0.5 pc, is lowered by a factor ~ 2 with respect 
to the more distributed population. An opposite trend 
would be expected in fragmentation models where the 
Jeans mass decreases with increasing gas density, leading 
to a corresponding increase in the abundance of BDs. This 
result seems to be best explained if a fraction of the dis- 
tributed population in Taurus is formed of low-mass stars 
and substellar objects ejected from the aggregates, either 
through dynamical encounters in a dense cluster (the colli- 
sion model) or through rapid dynamical decay in unstable 
sm all N-body systems (the em bryo-ejected model) . 

IKroupa fc Bouvierl l|2003a|) study the secular dynam- 
ical evolution of Taurus- Auriga type aggregates. They 
show that expansion occurs on timescales of a few Myr, 
corresponding to a few crossing times, due to two body 
encounters and gas removal. However, at the median age 
of the Taurus population (~ 3 Myr) the BD population 
does not separate significantly from the stellar popula- 
tion in these simulations. Long timescales are required, 
t > lOMyr, for a significant gradient in the relative 
number of BDs to develop within a central lpc sphere. 
Thus, secular dynamical evolution seems to require too 
long timescales to account for the pr esent spatial distri- 
butions of stars and BDs in Taurus. IKroupa fc Bouvierl 
( 2003b) have also investigated the case where BD ejections 
originate from dynamically unstable small N-systems (the 
embryo-ejection hypothesis) . The decay of small N groups 
occurs on time scale of a few 10 4 yr, ie the production of 
most BDs occur well before the end of gas removal, and 
the resulting escape velocities are larger than in the pre- 
vious scenario. The authors find th at the number of BDs 
per star seen in the Taurus survey of lBricefio" et al.N2002l) . 
taking into account BDs that escaped detection from their 
fields, and the central region of the ONC cluster are best 
explained if in both cases the embryo ejection model domi- 
nates and produces the same abundance of ejected BD per 
star of 25 % with a dispersion velocity < 2 km s _1 . The 
authors predict than in Taurus-Auriga, on scales ~ 1 pc 



from the center of the aggregates, there should be about 
twice as many unbound BDs than BDs bound to the ag- 
gregates. This scenario appears in good agreement with 
our observational results. 

Rece nt numerical simulations also seem to support this 
scenario. iGoodwin et alJ (2004a) follow the fragmentation 
of a distribution of 5 M© cores with turbulence properties 
similar to the Taurus cores. The resulting mass function 
is bimodal. Roughly 50 % of the sources, predominantly 
the low-mass ones, are ejected from the cores through 
rapid dynamical evolution before they could accrete signif- 
icantly, and form a flat low-mass distribution component. 
A high-mass log-normal component results from objects 
which remain bound in the core and grow by accretion. 
The ejected population comprises twice as many low-mass 
stars as substellar objects. Ejection velocities are 1-2 km 
s _1 almost independent of mass, leading to the same spa- 
tial distribution of BDs and low-mass stars. As mentioned 
before, in these simulations the abundance of BDs rela- 
tive to stars is on the order of 20 %, independent of the 
exact level of turbulence, and in close agreement with the 
Taurus observations. 

The line emission and infrared excess emission prop- 
erties of the very low mass population in Taurus indicate 
that accretion/ejection processes similar to the ones ob- 
served in the higher mass T Tauri stars proceed down to 
the substellar regime and that the inner disk lifetimes is 
similar in brown dwarfs and stars. These results have been 
taken by different authors as indication that brown dwarfs 
form through the same processes as stars and as an argu- 
ment against ejection models. However, the survival of the 
inner disk in young brown dwarfs may not be incompatible 
with the ejection model. Indeed with the typical inferred 
accretion rates of 10 -10 - lO -11 M0yr -1 , and with a typical 
disk mass of 1 % of the central stellar mass, survival times 
for an accretion disk trun cated at Rn„ t . = 10 AU around a 
50 Mj BD are a few Myr l|Liu et al.ll2003|) . ie on the same 
order as the age of the Taurus population. Viscous evolu- 
tionary time-scales appear to be on the same od er. Indeed, 
t u ~ ^(#) 2 (e.g. lArmitage fc Clarkel lll997ln scales as 
TV/- 1 / 2 w here M is the central stellar mass. Viscous time- 
scales for disks of R OM t = 10 AU around brown dwarfs of 
50 Mj are on the order of 2 x 10 6 yr (for a values of 10~ 3 
at r=10 AU). Thus accretion signatures may well last for 
a few Myr around brown dwarfs, even if small disks, with 
outer radii on the order of 10 AU, surround them, as ex- 
pected if they have been ejected early in their evolution 
from their birth site l|Bate et al.ll2002j) . 

In summary, the embryo-ejection model reproduces 
both qualitatively and quantitatively the spatial distribu- 
tion of stars and BDs in Taurus as well as the evolution of 
the 1Z SS ratio and does not appear incompatible with the 
accretion properties of the Taurus substellar population. 

6. Summary and conclusions 

We have presented in this work results from a spectro- 
scopic follow-up study of 79 very low mass Taurus can- 
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didates, selected from a large scale optical (I,z') survey 
covering a total area of ~ 30 square degrees. Our spec- 
troscopic survey is 90 % complete for photometric candi- 
dates later than spectral type M4V, 100 % complete for 
substellar candidates (spectral types later than M6.25V). 
Our corresponding mass completeness limit is 30 Mj for 
ages < 10 Myr and Ay < 4. We derive reddening values, 
spectral types and luminosity classes from a spectral fit- 
ting procedure and Na I equivalent widths measurements. 
We identify 17 new Taurus members among which 12 are 
brown dwarfs. We investigate their accretion properties, 
spatial distribution and abundance relative to stars. The 
main results are the following: 

1. A large fraction of the newly identified Taurus mem- 
bers show indication of accretion/outflow signatures. 
Two of the new VLM members (CFHT-Tau 19 and 
CFHT-Tau 21) display near- infrared excess emis- 
sion and optical forbidden lines indicative of accre- 
tion/ejection processes similar to the ones observed 
in the higher mass classical T Tauri stars. In addi- 
tion, 4 of the new substellar members (CFHT-Tau 6, 
CFHT-Tau 8, CFHT-Tau 11 and CFHT-Tau 12) ex- 
hibit Ha emission in excess of chromospheric levels 
and are likely CTTs substellar analogs. From levels 
of Ha emission, we estimate a fraction of accreting 
sources of 42 % in the Taurus substellar population, 
similar to what is observed in the higher mass T Tauri 
stars. 

2. We derive a new estimate of the substellar to stellar 



ratio in Taurus of 7?. s , = 



_ jV(O.O3-O.O8M ) _ 



JV(O.O8-1OM ) 



0.23 ± 



0.05, 1.3 times higher than the last estimate from 
iLuhmanl l)2004|) derived from spatially less extended 
studies. This value now appears fully consistent with 
the value derived with similar assumptions for the 
Tr apezium cluster pop ulation of TZ SS — 0.26 ± 0.04 
bv lBriceno et alJ l|2002l) and suggests an average pro- 
duction rate of one BD per 4 stars in both regions. 
3. We find a strong indication that the abundance of BDs 
relative to stars in the centers of the Taurus aggregates, 
on scales ~ 0.5 pc, is lowered by a factor ~ 2 with 
respect to the more distributed population. 

The similarity of the abundance of BDs with respect to 
stars in both the Taurus and ONC populations seem to fa- 
vor ejection models for the origin of the substellar popula- 
tion. The similar spatial distributions of the Taurus stellar 
and substellar populations but the detectable gradient in 
the BD abundance relative to stars at a characteristic age 
of ~ 3 Myr favor BD ejection through dynamical decay 
of unstable small N systems (the embryo-ejection model) . 
Further tests are however required to confirm this hypoth- 
esis, in particular multiplicity studies of the distributed 
population. In addition, our analysis relies solely on the 
relative abundance of BDs with respect to stars. Detailed 
and complete studies of the shape of both the substellar 
and stellar IMFs in star-forming regions of different envi- 
ronments are clearly required. 
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Fig. 12. The spatial distribution of the newly found Taurus BDs (sources wit h SpT < M6.25, triangle symbols) is 
shown superposed upon a Taurus CO map from llJngerechts fc Thaddeus (jl98lt (grey scale). Black dots indicate the 
distribution of previously known Taurus stars; filled white circles correspond to previously known BDs. The black 
curve s hows the outline of o ur o ptical CFHT sur vey, while the spatial coverage of the previous surveys from iLuhmanl 
( 200d). lBriceho et alJ ((20021) and lLuhmanl l|2004h is shown in dashed lines. 



